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Abstract 

The variation ranges of the axial B^o crystal-field parameters, for fc = 2, 4, 6, of tetragonal including 
cubic crystal-field Hamiltonians Ti-cp for all possible orientations of the relevant reference frame 
are studied. The distinguished 2-axis directions fixed by the maximal absolute values of B^q are 
analyzed. The diagrams for any tetragonal TicF parametrization depicting the maximal values 
of -^^^j where Mk is the 2'^-pole modulus, as a function of the x = or ratios, together 

with the distinguished directions are presented. The magnitudes and the relevant distin- 

guished directions are the discriminants of all the equivalent parametrizations. They vary within 
the intervals (0.7395, 1], (0.6074, 1] for tetragonal k — 4 and tetragonal k — 6 Hcf components, 
respectively. Such specified directions determine the mutual spatial orientation of the component 
2'^-poles of the Hcf, and due to their rigid coupling in the Hcf, they also refer to the global Hcf 
parametrization. This approach demonstrates the difference in fitting capability between the real 
and complex isomodular Hcf parametrizations. 
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1. Introduction 



In the paper by the equivalent TicF parametrizations we understand those referring to the same real 
crystal-field (CF) potential but expressed in variously orientated reference frames. Therefore, they can 
be reduced to the identical form through the appropriate rotations of the involved reference systems. 
However, the problem which generally remains, especially in the lower symmetry cases, is the lack of 
knowledge of the initial reference frame orientation in which the parametrization has been defined. 
The equivalent parametrizations can be identified and distinguished by their maximal axial parameters 
Bfco or their maximal absolute values, max|Sjto|, achieved in the specified reference frames. In general, 

max|i?;jo| ^ = \ Wkq\^ ( reach the values only under certain conditions for the crystal- 
field parameters (CFPs) ratios. It has been shown that the max ^^°^ magnitudes turn out to be 
useful discriminants of the equivalent parametrizations [1]. In the case of CF potentials of tetragonal 
and cubic symmetries expressed in their symmetry-adapted reference frames, there is no difficulty to 
establish such an initial frame orientation taking the four-fold axis as the ^;-axis. The only ambiguity 
resulting from the double choice of the x and y axes directions manifests itself solely in the change of 
sign of the tetragonal -Bfe±4 CFPs [2]. The absolute values of the |a;44| = |^| and |x64| = |^| ratios 
appear to be sufficient for the complete and explicit description of the parameterizations. 

There are in sum twelve point symmetry groups of the tetragonal and cubic symmetries. The 
parametrizations of the corresponding TicF^ split into two kinds [3,4,5]. The first concerns the tetrag- 
onal I [6] {D4, Civi D2d, D^ti) and the cubic (T, T/^, O, T^, Oh) point symmetries with purely real 
CFPs: 5 and 2, respectively. The second concerns the tetragonal II [6] (C4, S4, d^h) point symmetries 
which need nominally 7 CFPs, including the two pairs of complex ones. 

In compliance with the tetragonal point symmetry of the central ion its WcF parametrizations can 
include at most one pair of complex-conjugate CFPs per component multipole (except the quadrupole): 
-B44, -64-4; -B64) -B6-4- It means that the component multipoles parametrizations treated separately can 
always be reduced to a real three-parameter form. This reduction can be achieved by the appropriate 
reference system rotation about its z-axis [7]. 

Since the individual multipoles are characterized by the separate irreducible representations of the 
three-dimensional rotation group, the maximal B]^q variation ranges and the relevant distinguished di- 
rections are analyzed for each multipole individually. So, we use the real form of the T^cF parametriza- 
tions as the initial one. The more general complex parametrizations mentioned above are automatically 
included in the analytical approach based on continuous rotations of the reference frame. Formulating 
the global T^cf an additional problem can arise when all the three effective multipoles have to be ex- 
pressed simultaneously in a common reference frame (either known or a nominal one [8]). Then a pair 
of -Bjt4, -Bfc_4 CFPs (A; = 4 or 6) can occur in the tetragonal JicF parametrizations and an additional 
degree of freedom (the phase of the complex BkA CFP) is introduced into the fitting procedure. This 
problem is discussed in subsection 3.4 based on a classical example of ^4 symmetry. 

Rotating the initial reference frame by the (a,/?, 0) Euler angles within the ranges that cover 
the complete variation of the axial CFPs, two-dimensional functions Bko{a,P) are obtained. These 
graphs or their projections (contour maps) are characterized by the unique arrays of certain points 
further called the particular points, at which both the partial derivatives dBko/dp and dBko/da 
simultaneously equal zero [1]. These particular points can be classified as maxima, minima or saddle 
points by means of the higher derivatives either directly from the contour-line maps or from the 
three-dimensional BkQ{a,(3) graphs. It is recommended to present the maps in the normalized and 
dimensionless form dividing all their values by M^. According to the orientation of the parametrized 
CF potential with respect to the initial reference frame, the particular points arrange themselves into 
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the characteristic constellations. The values at these points can be degenerate. The maximum 

values, , are achieved in one or in several particular points. The plotted functions ^ 

(X44), {xq4) allow one to find directly the maximal values and the corresponding distinguished 

directions for any x ratios. Additionally, these practical schemes give also the values and directions 
for all the remaining particular points on the maps. It is worth emphasizing that the diagrams have a 
general character and make allowance for the explicit analysis of any tetragonal TLcf parametrization. 
Since the diagrams are of general character all the relevant dimensionless values are given with a fairly 
high accuracy of order of 10~^. 



2. General formulation 

Throughout the paper the tensor notation (by Wybourne [9]) for the T^cF parametrizations is consis- 
tently used. And so, in the case of the tetragonal I symmetries the parametrization has the form: 

7^CF(tetra I) = 520^^'^ + B^oCj,''^ + B44 {cf + C^_!l) + B^^cf^ + Bq^ {cf + ci']) , (1) 

whereas in the case of the tetragonal II symmetries the -6^.-1-4 are complex. In short, =J2 Cq{Oi, (pi) 

i 

where Cg{6i, ipi) are the normalized spherical harmonics [4,5,9] and the summation runs over i denoting 
all the open-shell central- ion electrons of the polar coordinates {9i,ipi). 

Here wc investigate the i?fco variation with the 2;-axis direction defined in an arbitrary reference 
frame by only two Euler angles (a, f3, 0) with respect to the initial frame [10,11], where the third angle 
7 of the rotation about the final z-axis may be assumed to be zero as an arbitrary in the case. The 
transformed Bj^q are related to the original B^q as [1,5,10,11] 

k k 

B'ko = E <^(«,/3,0)i?fe, = E (/3,a)5fe„ (2) 

q=—k q=—k 

where vll;^{a,P,0) = cf^(/3, a) arc the middle row matrix elements of the P'^'^'(a, /3, 0) rotation 
matrices [10,11]. Below in the paper B^q correspond to the original reference frame, whereas their 
primed counterparts stand for the rotationally transformed ones. 

Naturally, all the equivalent TCcf parametrizations have to be characterized by the same variation 
ranges of their CFPs including the axial ones, in other words by the ranges with the same upper and 
lower limits. These particular maximal and minimal values can serve as useful discriminants for the 
whole classes of the equivalent parametrizations. In order to standardize the equivalent parametriza- 
tions the magnitudes can be recommended as the best discriminants. Consequently, the 
direction of the transformed reference frame z-axis relevant to these discriminants are now the com- 
mon distinguished directions for all the equivalent parametrizations. This approach is especially useful 
in the case of low-symmetry CF potentials [1]. 

The axial CFPs play an extraordinary role since they are real and invariant relative to rotations 
about the z-axis. They are also independent of the z-axis inversion, whereas the remaining CFPs 
become then inevitably transformed according to the rule Bj^q — > Bi^_q, as the effect of the rotation 
by (0,7r,0). 

The relationship between 5jJ.q CFPs and the two rotation angles {a, (3, 0) of the reference frame can 
be expressed in the form of two-dimensional contour maps "^''^""'^^ with iso-^^ lines joining the points 

of identical -j^ values. For the equivalent parametrizations such maps can be always transformed into 
each other. The whole ranges of the rotation angles variation < a < 27r and < /3 < tt can be 
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limited according to the TicF symmetry. For the tetragonal symmetries (I and II) with the z-ajds 
along the four-fold axis the maps are periodical functions of a with the period of tt/2. Furthermore, 
since -B^q CFPs are the same at points {a, (3) and (a + vr, vr — what results from their independence 
of the z-axis orientation, the range of /3 can be confined to the interval (0, 7r/2). For the clarity reasons 
the maps are presented within the ranges < a < tt and < /? < tt. 

As the paper concerns exclusively the rotational transformations of the axial -B^q CFPs, only the 
middle rows of the relevant rotation matrices T>^^\a,l3,G) [10,11] are employed. All the non-axial 
CFPs transform themselves according to the remaining rows [5,10,11]. After the initial reference 
frame rotation the T^cF form obeys the central-ion point symmetry with respect to the new frame. 
Additionally, in the distinguished reference frames some correlations between the CFPs can occur, like 
that found by Burdick and Rcid [12]. They proved analytically that the -B^ii CFPs become zero at 
the maximum, minimum and saddle points of the Bi;o{a,f5) functions. 



3. The axial components in equivalent parametrizations of tetragonal 
and cubic Hcfs as functions of ^ or ^ ratios 

-D40 -D60 

3.1 The 2^-pole ctxial parameter B20 in cirbitrarily orientated reference frame 



In cubic potentials the quadrupolar moment is completely compensated and, in consequence, its 
components do not appear in any reference frame. Hence, in the most general case the cubic T^cF can 
be parametrized by at most 22 instead of 27 CFPs. 

For the tetragonal potential expressed in the symmetry-adapted system there exists only the axial 
term B2qCq^ . In any other reference frame rotated by a and (3 one gets [1,10,11]: 

B'2^ = S2oPg^(a,/3,0) = {-iycf\^,a)B2o, (3) 

B' 

and the ratios for different q are correlated. For ^ = 

^20 = \ (3c' - l) i^20, (4) 

where c = cos/? and thus i?2o depends solely on (5. The derivative dB2Q/df3 = —3csB2o, where 
s = sin/3, and is equal to zero for /3 = 0,7r or tt/2. For these /? values i?20 J^eaches its extreme 
magnitudes: B20 and — ^B2o, respectively. Hence, the distinguished direction for which '-j^ achieves 
1 is the symmetry axis, whereas in the symmetry adapted system equals 1 automatically. For 
P = arccos which is the angle between the edge and diagonal in a cube, B2Q vanishes, but then 

the four remaining CFPs appear: I-B22I = |-B2_2| = "^-620 and I-B21I = |-B2_i| = ^^20. 

3.2 The 2^-pole axial parameter B40 in arbitrarily orientated reference frame 

The fourth-order TCcF term in the tetragonal I and cubic CF potentials described in the symmetry- 
adapted system has the form [3-5,13] 

nS^ = B,o[c^'^+x,{ci'^+C^^l)], (5) 

and the modulus of the 2^-pole is equal to 

M4 = |54o|(l + 2x1)1/2. (6) 
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For the cubic case 3^4 = = 0.5976 [6,7,13]. Eq.(3) yields for any rotation by a and /3 the relation: 
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3 - SOc^ + 35c^ + Xi VtQs^ cos 4al B^q . (7) 



Therefore, 

5^0 3- 30c2 + 35c'^ + X4V70s^cos4q! B40 



M4 8(1 + 2x1)1/2 1^40 r 

Using Eq.(9) for any 0:4 the map -^(a,/?) can be depicted within the required ranges of a and (5 
angles. The partial derivatives 



1 dB'^Q _ cs(15-35c2 +X4\/70s2cos4a) B. 
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M4 dp ~ 2(1 + 2x1)1/2 IS40I 

1 55^0 _ _a;4\/70s^sin4a B^q 
M~4~d^ ~ ~ 2(1 + 2x1)1/2 ■ \B^\ 

simultaneously vanish in the following five cases: 



(9) 



a) s = 0, 

b) c = 0, cos 4a = 1, 

c) c = 0, cos4a = — 1, 

d) 15 - 35c^ + X4\/70s^ = 0, cos Aa = 1, 

e) 15 - 35c2 + X4^s^ = 0, cos 4a = -1, (10) 

IB' 

which correspond to the following points on the map -^(a,/3) 

a) P = 0, 

b) p = —, a = 2n— , 

c) /?=|, a = (2n + l)^, 

d) /5 = arccosf^l^^V'^ " = ^4 
^ V^4V70 + 35; 4 

, ^ /'x4\/70 — 15\ , .TT , , 

e) p = arccos = , a = (2n + l) — . (11) 

^ V^4V70-35y ' ^ M ^ ^ 

In Eqs (11) n stands for any integer. Let us denote by ji^jg |, i = a . . .e, the absolute value of B^q for 

the {i)-th condition in Eq.(lO). Since the plots -^^{x^) for (a) and for the pairs (6, c) and {d,e) are 
symmetrical with respect to X4 (Eqs 8, 11), only the positive values of X4 are taken into account. 

For each condition (Eqs 10, 11) the expressions -j^{xa) can be represented by the following plots 
(Eq. 8): 

a) ^(X4) = (1 + 2x1)-^/^ 

For X4 = 0, i.e. for the purely axial T^CF, this function reaches maximum 1. For the cubic case 
(x4 = 0.5976) it amounts to 0.7638, and monotonically falls to zero when X4 — 00 (Fig. 2). 
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For 0:4 = this function equals to 0.3750. For the coplanar square ratio (X4 = 1.3944) it reaches 
its maximum 0.8292, and then asymptoticahy decreases to 0.7395 when X4 — ^ 00. 

This function for X4 = takes the same value as in (b) since they are symmetrical. When X4 
rises it drops to zero at X4 = 0.3586, and then monotonically increases to 0.7395 for X4 — > 00. 

d) The condition < < 1, where = is fulfilled for all X4 > because X4 > -1.7928 

is sufficient. Substituting the above expression for and cos 4a = 1 into Eq.8 yields 

^ _ 1^4770-151 
^4 ^^^^ - (,,y70+35)(l+2.|)'/^ 
For .T4 = the function equals 0.4286, it has zero-point at X4 = 1.7928 and then it pursues zero 
for ,T4 00 passing through the weak and flat maximum 0.0489 at X4 = 5.1532. 

e) The condition < < 1, where = f^^^rH' confines its domain of determinancy to 
X4 < 1.7928. Within this range 

- X-4V70+15 
^4 ^"'"^ " (|.4x/70-35|)(l+2.1)'^^- 

For X4 = this function equals 0.4286, and for the limiting value X4 = 1.7928 it is 0.5504, i.e. 
the same as in (c) at that point. Within the interval < X4 < 1.7928 the plot passes through 
the subtle local maximum of 0.5092 at X4 = 0.5977, and minimum 0.5065 at X4 = 0.9701. 

The plots i^(x4), i = a,b,c,d and e corresponding to the five conditions (Eqs 10, 11) are presented 

in Fig. 2. The diagram shows the positions and values of all particular points on the 3^(ct,/3) map 
for any tetragonal HcF as a function of X4 (see Figs 5-7) . Any vertical straight line drawn for a given 
X4 intersects the diagram branches (from a to e) in three to five points (depending on the X4) directly 
determining the distinguished directions and values including the 

As is seen from the diagram the four- fold axis, i.e. the initial z-axis direction ([001]), is the 
distinguished one for |x4| < 0.5976. For all the remaining X4, i.e. |x4| > 0.5976, this distinguished 
direction lies along the square diagonal [110]. The crossover point corresponds exactly to the cubic X4 

value. The ""^f^"' value varies within the range from 1.0000 (for X4 = 0) down to 0.7395 for X4 00, 
achieving on the way a local minimum 0.7638 for the cubic X4 value, and a local maximum 0.8292 for 
the coplanar square X4 value. 

3.3 The 2^-pole cixial parameter Bqq in cirbitrarily orientated reference frame. 

In turn, the sixth-order T^cF term in the tetragonal I and cubic CF potentials in the symmetry adapted 
system has the form [3-5,13]: 

4V^ = 56o[C7r+xa(cf (12) 
and the modulus of the 2^-pole is equal to 

Me = |S6o|(l + 2x^)1/2. (13) 

For the cubic case xq = = 1.8708 [6,7,13]. In the reference frame rotated by (a, (3,0) the axial 
CFP Bqq transforms to (see Eq.3): 

B'^Q = -L [_5 + lOSc^ - 315c^ + 231c^ + x^^VUs'^i-l + llc^) cos 4a] Beo- (14) 
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Hence, 



B'qo -5 + 105c2-315c^ + 231c6 + X63^s^(-l + llc2)cos4a 



-'60 

Me 



60 



16(l + 2xi)V2 



5, 



60 



(15) 



The particular points of the map ■j^{a,/3), especially the point (or points) corresponding to the 
'"^f^"^ , are those for which both the partial derivatives dE'e^/dp and 83'^^/ da simultaneously become 



zero 



1 dBL 



60 



Me dp 



3 -35 + 210c2 - 231c^ + X6Vl4(-13 + 460^ - 330^) cos 4a Bqq 



1 SB' 



60 



(1 + 2xi)V2 

3 X6yi4s^(-1 + llc^) sin4Q; Bqq 
4 



\B, 



60 1 



Me da 4 (l + 2x2)V2 iSeol' 

There are six possibilities of simultaneous vanishing of the partial derivatives 
a) s = 0, 

h) c = 0, cos 4a = 1, 
c) c = 0, cos 4a = — 1, 



(16) 



cos 4a 



0.5238 



d) c2 = 

e) (231 + 2,'iy/lAxe)c^ - (210 + 46Vl4x6)c^ + (35 + 13\/l4x6) = 0, cos 4a = 1, 
/) (231 - 33\/l4xe)c'^ - (210 - 46^X6)0^ + (35 - 13\/l4xe) = 0, cos 4a 
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.4 /oin I /i~T^ _ ^ „2 



1. 



(17) 



The last two conditions, as the functions of Xe yield the double solutions, i.e. four solutions for 
xe > 0.7195 and three ones for |xe| < 0.7195, so in sum there are eight or seven particular points 
on the map apart from their symmetrical images due to the presence of the four-fold axis as well the 
invariance of the axial CFPs with respect to the z-axis sense. 

The conditions specified above in Eq.l7 (from a to /) correspond to the following particular points 
on the map ^(a,/3): 



TT 

^4' 

n + 1)-, 



a) 


P = 


0, 


b) 




TT 

2' " 


c) 


P = 


vr 

2' " 


d) 


P = 


72.45°, 


e) 


13 = 


arccos 


e') 


P = 


arccos 


/) 


f3 = 


arccos 



1 

4 



-0.5238\ 
X6 J 



/') P = arccos 



1 


15^ 


/li + 46x6 - 2(100xi + A8^ 


/lixe 


+ 210)^/2 


33 










1 


15\ 


/14 + 46x6 + 2(100xi + A8^ 


/lixe 


+ 210)1/2 


33 




\/T4 + 2x6 






1 


15\ 


/li - 46x6 - 2(100xi - 48\ 


/lixe 


+ 210)1/2 


33 




\/l4-2xe 






1 


15\ 


/14 - 46x6 + 2(100x^ - 48\ 


/lixe 


+ 210)1/2 



33 



14 - 2x6 



1/2 



1/2 



1/2 



1/2 



„ TT 

a = 2n— , 



a = 2n 



TT 

4' 



a 



a = 



(2n+l)|, 



TT 



(2n+l)-. (18) 
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In the cases (d), (e) and (/) the definite domains of determinancy of the involved expressions have to 
be taken into account. In consequence, the xq ratio has to fulfil the following conditions: \xq\ > 0.5238, 
xq > —0.7195 and xq < 0.7195, for (d), (e) and (/), respectively. 

Unlike for the particular points in (a — c) cases their location in (d — f) cases does depend on the 
xq ratio. The (d) case turns to be quite intriguing. The corresponding particular points (Eq.lSd) lie 

along the straight line /? = 72.45°, and the -j^{xe) function takes along it the constant value ^j^^^i-^a 
regardless of the a rotation angle. 

In the (e) case the /3 angle varies from 90° for the limiting value xq = —0.7195, through 62.04° 
for xq = 0, and monotonically goes down to 51.12° for xq oo. The P{xq) plot for the (/) case 
is symmetrical to the previous one with respect to xq. Similarly, the symmetrical plots for the (e'), 
Eq.(18e') and (/'), Eq.(18/'), cases are characterized by /? = 51.12° for xq Too, /3 = 33.88° for 

Xq = 0, and /? = for — d=oo, respectively (Fig.3). Since the plots {xq) for (a), (d) and the 
pairs (6, c), (e, /) and (e', /') are symmetrical with respect to xq (Eqs 15, 18) only the positive values 
of xq ratio are considered. 

For the conditions defined above (Eqs 17,18) the expressions -j^{xq) reduce themselves (Eq.l5) 
to the following plots: 

a) ^(x6) = (l + 2xi)-V2 

For xq = the maximal possible value 1 is reached, for the cubic ratio (xq = 1.8708) it amounts 
to 0.3536, and then monotonically falls down to zero as xq ^oo. 

L>j Me y-^^J - 16(l+2a:2)l/2 

For Xq = the function is equal to 0.3125, for the coplanar square ratio {xq = 1.1225) it reaches 
its maximum 0.5863, and then slowly falls to the limit 0.4961 as xq — oo. 

r\ \^(Ta) - |-5 + 3v^X6| 
^) Me y-^<i) - 16(l+2a;2)i/2 

As above, for xq = the function equals 0.3125, and when xq rises it falls to zero for xq = 0.4454, 
then monotonically increases to the asymptotic value 0.4961 if ^ oo. 

A\ 1^6? I (^^\ - 0.1322 

Me V-^6J - (1+23:^)1/2 

Here, the domain is limited to \xq\ > 0.5238, (Eq.l8). For xq = 0.5238 the function amounts to 
0.1062 and monotonically tends to zero for xq oo. 

X ISgoii/' \ _ |-5+105ti(a;6)-315t2(a;6)+23U3(a;6)+X637l4(l-ti(3:6))2(-l+lUi(a:6))| 
^) Me \^(^) - 16(l+2a;2)i/2 ' 



, . / \ 1 15v^+46a;6-2(100x2+48v^a:6+210)i/2 , nT-mc 

where ti{xe) = — - — ^_^2xe — — ' ^ — 

For Xq = the function equals 0.3321. For the cubic ratio {xq = 1.8708) it reaches quite a subtle 
maximum 0.6285, and then slowly falls to 0.6074 if xe — oo. 

IB^o'l f x _ \-5+105t2{x6)-3154ix6)+23iq{xe)+xe3Vu{l-t2{xe))'^{-l+llt2{xe))\ 



e ; Me - 16(1+212)1/2 

Where t^ixe) = i i5x/i4+46.a+2(ioo.i+48x/n.e+2io)V2 



33 Vu+2xe 

For Xq = the expression equals 0.4148, it has zero-point at xq = 1.2481, next slightly rises to 
the maximum 0.0710 at xq = 3.6498, and then slowly falls to zero if xq oo. 

fN / X _ |-5+105t'i(a;6)-315t;2 (xs)+231tf (x6)-a:63\/l4(l-t;(a:6))2(-l+llf;(x6))| 

^) Me Va^ej - 16(1+2x2)1/2 

Where t[{xQ) = l 15VT4-46.a-2(1^.2 48VT4.e+210)V2 ^ ^ ^^^^ 

^ ' "^o y 14:— 2XQ 
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For xq = the function equals 0.3321. At the hmiting point {xq = 0.7195) it amounts to 0.1348, 
and within the interval (0 < xe < 0.7195) the plot passes through the minimum 0.1050 at 
Xd = 0.5595. 

f/A l^go^lc ^ _ |-5+105t^(x6)-315t^^(3:6)+231f^3(a;6)--a:63\/l4(l-f^(a;6))^(-l+llt^(x6))| 
I ) Mfi y^^l - 16(1+21^)1/2 ' 

wVi^r^ +'(^ \- 1 15%/n-46rr6+2(100a;g-48N/nx6+210)V^ 

wnere t2[xe) - ^ Vu-2x6 • 

For Xq = the function equals 0.4148 and then monotonically rises to 0.6074 as — oo. 

The j^" (xg) plots, where i = a,b, . . . , f , for all the eight conditions (Eqs 17,18) arc presented in 
Fig. 4. The straight line corresponding to any chosen xq ratio intersects the diagram branches in four 
to eight points depending on the xe value (see Figs 5-7). The points of intersection determine the 

positions of the particular points and the respective values including the —j^^{xq) discriminant. 
The characteristic degeneracy of some particular point values in the diagrams manifesting itself as the 
intersection points of their different branches (Figs 2,4) reflects some definite symmetry encoded in 
the xq ratio. 

As is seen from the diagram (Fig.4) the four-fold axis, i.e. the initial z-axis direction ([001]), is the 
distinguished direction for \xq\ < 0.9179. For all the remaining xq, i.e. \xq\ > 0.9171, the distinguished 
directions (roughly [111]) and the corresponding maximal values are determined by the e-th branch 
of the diagram. In the cross-over point (xq = 0.9179) this maximal value amounts to 0.6103. For 
Xq > 0.9179 the distinguished direction forms with the reference frame z-axis the angle /? given by 
Eq.lSe. In the crossover point this angle amounts to 56.69°, for the cubic ratio (xg = 1.8708) it equals 
54.74° (exactly [111] direction), and then asymptotically pursues the limit 51.12" as xg oo. 

The —^^^{xe) can vary within the range from 1.0000 (for xq = 0) to 0.6074 (for xg —>■ oo) 
passing on the way through the local minimum 0.6103 (for xg = 0.9179) and the maximum 0.6285 (for 
Xe = 1.8708). It is worth noticing that the distinguished direction crossover takes place for xg = 0.9179 
which is smaller than both the coplanar square ratio 1.1225 and the cubic ratio 1.8708. 



3.4 Tetragonal II Hcf with complex parameters — the maiximal axial CFPs in arbi- 
trarily orientated reference frame 

The tetragonal II WcFS containing two pairs of 5^4, -Bfc-4 CFPs (fc = 4, 6) in the symmetry adapted 
system deserves more careful analysis. If only one pair of the complex CFPs occurs in the T^CF, there 
always exists the equivalent real parametrization. The complete parametrization of the tetragonal II 
TicF contains seven instead of five CFPs. Taking into account the rotational properties of TCcf it 
is recommended to use the pairs of complex-conjugate parameters: = |-BA;4|e^^'^'=* and -Bfc-4 = 
|-Bfe4|e~'^'^*4 rather than the equivalent pairs ReSfe4 and ImBk^. Between the former and latter pairs 
of the CFPs the following relations hold: Bk4 = KeBk^ + ilmBk^, Bk-4 = KeBk4 — ilm5fc4, \Bk4\ = 
[(Re-Bfc4)^ + (Ini-Bjfc4)^] (fk4 = \ arctan {["^''"^ [4,7]. Rotation about the z-axis (the four-fold axis) 
by the angle ^pkA eliminates the lmBi.4. Then, the new real B'j^^ is equal to \Bf;4\, and the tetragonal 
parameters for the remaining k transform correspondingly, whereas the axial CFPs (g = 0) remain 
unchanged. Therefore, the tetragonal II T^cF for d-electron transition ions can always be parametrized 
purely by the real CFPs. Similarly, for the hexagonal C^h symmetry the number of CFPs can be 
reduced from five down to four and is equal to that for D^^ symmetry, when the Imi?66 is absent due 
to the higher symmetry [7,8,14]. Such elimination of one of the two lm.Bk4 CFPs {k = A or 6) fixes 
the reference frame, and is equivalent to the lost of one degree of freedom, i.e. the rotation of the 
reference frame [6]. This fixed frame is imposed on the remaining effective multipoles in the T^cF to 
be parametrized. 
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Let us consider as an example the T^cF of ^4 point symmetry occurring in paramagnetic ions in 
the scheeUte type crystal structure [7,15]. Then, the complete seven-parameter HcF can be written 
in the form 

ncF{S4) = B^oCf^ + 54oCi^^ + |B44|e^^'^«C7f + |B44|e-^^'^«ci'i + 

+ B6oCi^^ + |S64|e^^'^64c{6) + |5g^|e-i4<^64c(6^^ ^^g^ 

where the Bj^^ and i?fe-4 CFPs are only conjugate (Eq.(l)). 

This form is usually reduced to the six-parameter one. The ^4-4 is eliminated rotating the reference 
frame by the angle ipn, and so we get 

HcFiSA) = B-.oCf'^ + B^qC^^^ + \Bu\{cf + C^fi)+B^oCf^ + 

+ |564|e*^('^64-cp44)(-_(6)^|^^^|g-i4(c^64-C^44^(6)^ ^20) 

This formula differs from the real one of the same M in one additional parameter (Eq.(l)) - the phase 
difference ((/?64 — ¥'44) of the complex CFPs, which is an extra parameter in fitting procedures. Now 
the max 1 5go I magnitude for all possible orientations of the reference frame is the maximum of the 
expression (Eq.(14)) 



\B' 



601 



1 

16 



-5 + lOSc^ - 315c'' + 231c^ + XQiiVlis\-l + llc^) cos 4 [a - {ipQ^ - ip^^)] ■ \Bqq\, (21 



where xqa = Thus the shift (9954 — ^aa) does not change the max|i?gQ|, but influences the relevant 

distinguished direction compared to that for the real parametrization and, in consequence, modifies 
the spatial orientation of the considered multipole. 

Let us note that in the case of the tetragonal I T^CF, i-e. for real CFPs, there is no such possibility 
- the orientation of the considered 2*^-pole is absolutely determined by the Xfc4-ratio. The result is 
that the complex HcF parametrizations enable some extra modifications of the spatial orientation of 
the component multipoles in fitting procedures when the absolute values of the CFPs are conserved. 



5. Discussion 



The maximal admissible magnitudes of the B^q CFPs, their variation ranges as well as the relevant 
distinguished directions depending on the reference frame orientation for all classes of the equivalent 
parametrizations for the tetragonal including cubic Hcfs have been thoroughly analyzed. The anal- 
ysis was inspired by observation that in any reference frame the entire magnitudes of cannot be 
generally achieved by the parameters of the TicY 2'^-pole component, particularly by the axial param- 
eters B^q for even k > 2. The only exception is the purely axial TCcps for which — = 1 if the 
reference system z-axis is in line with the symmetry axis. The quadrupole {k = 2) is a special case 
(not to mention vectors with A; = 1) since |i?2ol can reach M2 for X2 = -B22/-B20 = ^ [16], but this 
case does not concern the considered symmetries. 

The calculated dependencies of the — 1^*^(2;^) for k = 4, 6 on the TicF 2'^-pole composition 

represented by the appropriate x ratio are shown in Figs 2 and 4. They offer the general review 

B' 

of all particular points of the maps -j^{a,P) establishing their positions and giving the relevant 

values for any X4 and xq ratios. The values can serve as discriminants of all the equivalent 

parametrizations which have the same form in the common symmetry adapted system. 

The degree of reduction of the discriminants with respect to 1 is more tangible for A; = 6 

than for A; = 4. The whole ranges of the discriminants as functions of x are limited. The lower limits 
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of "^"1^^"^° ^ and ^ amount to 0.7395 and 0.6074, respectively, whereas for comparison in the case 

of quadrupole ^-^^ > 0.8660 [16]. The envelopes vs. x are characterized by a sharp local 

minimum at the intersection point of the two plots representing the two various conditions for the 
extrcmum, and a weak flat maximum for the ratios 3:44 = 1.3944 and xqa = 1.8708 corresponding to the 
co-planar square and cubic coordinations, respectively. These particular x ratios can be linked up with 
the coordination geometry. Taking into account the point symmetry of the central ion surroundings 
and 1/R dependence of the individual ligand potential one gets [4] 

-B44 ^ Vjo 

B40 2(2 cot^ e-12 cot2 + 3) ' 

Bu ^ ^^VU 5cot^g-l 

Beo 2 2cot6 0-3Ocot4 + 45cot2 0-5)' ^ ' 

where 6 and R stand for the spherical coordinates of the ligands forming a square. The sign of the 
ratios depends on the x and y axes orientation with respect to the square. 

I-B' I 

The distinguished directions corresponding to the maximal magnitudes of (x) are presented in 
Fig.8. All the directions are defined with respect to the reference frame in which the initial parametriza- 
tion is expressed [1,8]. 

The above investigations provide the distinguished directions for all the component multipolcs 
separately, but they are also obligatory for the resultant Hcf- These are the directions that could 
be used in a possible reduction of the number of independent CFPs [12], as well as examining the 
standardization of equivalent parametrizations [17,18]. The recognition of the distinguished directions 
for individual 2'^-poles gives their autonomous spatial orientation and throws light on the mechanism 
how the component multipoles combine into the resultant effect. The role of the mutual orientation of 
the component multipoles is clearly noticeable while the real (five-parameter) tetragonal parametriza- 
tions of Hcf are compared with the more general complex (seven- or six-parameter) iso-modular 
parametrizations. The latter is more capable in fitting procedures what implies directly from the 
possibility of modification of the component multipoles mutual orientation. Dealing with the lower 
symmetry Hcf, in which the component multipoles contain more than one complex-conjugate pairs 
of CFPs, some intra-multipoles modifications are also possible. 

The -j^{a,(5) and -^(a, /3) maps for three representative, well-documented parametrizations of 
the tetragonal Hcf^ significantly differing in rc44 and rcQ4 ratios are shown in Figs 5-7 as an example. 
The maps in Fig. 5 show the variation of the B^q and Bqq CFPs depending on the reference frame 
z-axis direction for the Pr^+ ion embedded in the KY3F10 matrix [19]. Its HcF of C^y point symmetry 
is characterized by the following parameters (in cm~^): B20 = —589, B40 = —1711, -B44 = 519, 
Beo = 657, Bqa = -196; M2 = 589, M4 = 1862, Me = 713; X4 = -0.303, xq = -0.298. 

The maps in Fig.6, in turn, apply to the cubic Hcf with the characteristic ratios X4 = 0.5976, 
and xq = —1.8708 [3-5,13]. The system Pr^+:Cs2NaRCl6 is given here as an example, and, the Oh CF 
potential is defined by the CFPs (in cm"^): B40 = 2326, B^o = 247; M4 = 3045, Me = 699 [20]. 

Successively, the maps in Fig. 7 refer to the B4Q and Bqq CFPs variation for Pr"*"^ ion in Pr2Cu04 
[21], i.e. in the €4^ CF potential of the ratios X4 = —0.830, xq = 9.99, and the CFPs (in cm^^) are 
correspondingly equal to: B20 = -294, B40 = -2302, B44 = 1910, Bqo = 147, B^ = 1469; M2 = 294, 
M4 = 2991, Mq = 1476. With the exception of the cubic case the remaining X4 and xe ratios come 
from fittings experimental data. 

The particular points marked on the maps (Figs 5-7) by the letters are also noticeable on the 
diagrams (Figs 2 and 4) as the intersection points of the appropriate diagram branches with the 
vertical straight lines corresponding to the involved X4 and xq ratios. The intersection points with the 
highest branches refer to the "^'^^^'^^ discriminants. 
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Figs 2 and 4 present the variation of the particular points location on the maps -^(«, P) with X4 

and xe- Since the ivlf P) maps can be considered as the representations of the parent parametriza- 
tions, further conclusions concerning the links and genealogy of various CF potentials, their optimal 
higher-symmetry approximations, etc. may be drawn based on the evolution of such maps while 
descending in the central ion symmetry, what can be simulated by its shifts inside the coordination 
polyhedron in various directions. Although the results presented in this paper concern the equivalent 
HcF parametrizations, from the point of view of the maximal axial parameters of component multi- 
poles and their relevant distinguished directions they have more general character going beyond the 
crystal-field approach. 
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FIGURE CAPTIONS 

FIG.l The dependence of on the reference frame rotation angle /? for tetragonal T^cfs (Eq.4). 

I^(»)| 

FIG. 2 vs \xi\ plots, where i = a,b,c,d,e denotes one out of the five conditions (Eq.lO) or 

the particular points of the map (Eq.ll). The vertical straight lines refer to: \x4\ = 0.303 

for Pr3+:KY3Fio [19]; \x4\ = 0.598 for Pr3+:Cs2NaRCl6 [20], and |a;4| = 0.830 for Pr2Cu04 [21], 

respectively. 

FIG. 3 P{xe) functions for (e) and (e') conditions (Eq.18); for (/) and (/') conditions the plots are 
symmetrical with respect to xq. 

FIG.4 '-^^ vs \xq\ plots, where i = a,b, c, d, e, e' ,f, f denotes one out of the eight conditions (Eq.l7) 

or the particular points of the map — (Eq.18). The vertical straight lines refer to: l^el = 0.298 
for Pr3+:KY3Fio [19] and \xq\ = 1.871 for Pr3+:Cs2NaRCl6 [20]; for Pr2Cu04 xq = 9.99 [21]. 
FIGS 5a, 5b The maps ^^^^^ and ^^^^ for Pr3+:KY3Fio [19], point symmetry; the initial 
CFPs (in cm-i) are B40 = -1711, B^o = 657; M4 = 1862, Mg = 713; X4 = -0.303, xq = -0.298; the 

particular points arc marked by the letters. 

FIGS 6a, 6b The maps and ^'""ms^^ for Pr3+:Cs2NaRCl6 [20], Oh point symmetry; the initial 

CFPs (in cm'^) are B40 = 2326, Bqq = 247; M4 = 3045, Me = 699; X4 = 0.5976, xq = -1.8708; the 
particular points are marked by the letters. 

FIGS 7a, 7b The maps ^^^^ and for Pt^+ in Pr2Cu04 [21], 64^ point symmetry; the 

initial CFPs (in cm-^) are B40 = -2302, Bqo = 147; M4 = 2991, Mq = 1476; X4 = -0.830, xq = 9.99; 
the particular points are marked by the letters. 

FIG. 8 The ranges of the distinguished directions characterized by the maximal values of jS^ol CFPs 
depending on the |x| ratios for tetragonal TCcp's (|x4| = 0.5976 and \xe\ = 1.8708 are the cubic ratios, 
the /? angle between the distinguished direction and the z-axis is given within the round brackets). 
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Fig.3. 
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Fig. 6a 
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Fig.7b 
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